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Abstract 

Nucleosynthesis calculations require nuclear level densities for hundreds or 
even thousands of nuclides. Ideally one would like to constrain these level 
densities by microscopically motivated yet computationally cheap models. A 
statistical approach suggests that low moments of the Hamiltonian might be 
sufficient. Recently Zuker proposed a simple combinatoric formula for level 
densities based upon the binomial distribution. We rigorously test the bi- 
nomial formula against full scale shell-model diagonalization calculations for 
selected sd- and p/-shell nuclides and also against Monte Carlo path integra- 
tion calculations. We find that the fourth moment is as important as the third 
moment to a good description of the level density, as well as partitioning of 
the model space into subspaces. 
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I. INTRODUCTION 



Nuclear level densities are important for the theoretical estimates of nuclear reaction 
rates in nucleosynthesis. The neutron- capture cross sections are approximately proportional 
to the corresponding level densities around the nuclear resonance region. The competition 
between neutron-capture and (3 decay determines the fate of the s- and r-processes. The 
abundance of s-process nuclei with nonmagic neutron number is inversely proportional to 
the neutron-capture cross section i.e., proportional to the level density. The waiting 
points of the r-process are determined by the balance between the rates of neutron-capture 
and of photoejection of a neutron. 

Reliable estimates of nuclear abundances require accurate level densities. Authors like 
Gilbert and Cameron gave uncertainties beyond a factor of 10 in their early calculations 
and uncertainties of about a factor of 8 and 5 are present in calculations of § and [Q, 
respectively. What is desired is a reliability within a factor of two. 

At low excitation energies (up to 5 MeV), level densities are usually extracted experimen- 
tally via direct counting from neutron/proton resonance data 0. For intermediate energies 
(between 5 and 15 MeV) they are extracted from charged particles spectra while for 
still higher energies the level densities are determined using the Ericson fluctuation analysis 
0. Recently, the Oslo group has reported on a new method to extract level density from 
primary 7-ray spectra without the assumption of any model P]. 

Such experimental measurements are very labor intensive, and reaction network calcu- 
lations require hundreds or thousands of cross sections, sometimes for unstable nuclides. 
Therefore one turns to theoretical models for help. 

Most conventional calculations of the nuclear level density are based on the Fermi gas 
model. The most widely used description of the nuclear level density is the Bethe formula, 
based on a gas of free nucleon 0. This parameterization works well for the level densities 
of many nuclides at low energies and in the neutron energy region [|10| especially when 



one uses the modified "backshifted Bethe formula" [11,12]. Unfortunately the fitted single- 
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particle level density parameter and the backshift parameter are nucleus-dependent and are 
not derived theoretically. Rauscher and collaborators presented a global parametrization of 
nuclear level densities within the back-shifted Fermi-gas formalism |T^. They employed an 
energy-dependent level density parameter and a backshift parameter to attain an improved 
fit of level densities at neutron-separation energies. Despite all this, the fact remains that 
the input parameters are just that: empirically fitted parameters. 

Can one look to theory for help? The interacting shell model and other microscopic mod- 
els accurately describe spectra and transition for a broad range of nuclides. On the other 
hand, "traditional" shell-model codes diagonalize the Hamiltonian in a large-dimensioned 
basis of occupation-state wavefunctions but use the Lanczos algorithm to extract only a 
handful low-lying states. The level density requires all eigenstates and thus complete diag- 
onalization, a computationally forbidding requirement. 

An alternative to diagonalization is the Monte Carlo path-integral technique [|1^], which 



is well suited to thermal observables |T5HT7[]. Nakada and Alhassid extracted both the 



single-particle level density parameter and the backshift parameter from the microscopic 
Monte Carlo densities by fitting them to the backshifted Bethe formula and also found 
some shell effects in their systematics. Recent developments include the particle number 
reprojection method flSl, which can calculate level densities of odd-even and odd-odd nuclei 



despite a sign problem introduced by the projection of an odd number of particles. Path- 
integral methods are limited, however, to interactions that are free of the 'sign problem.' 
Therefore we consider alternatives. 



One possibility is nuclear statistical spectroscopy [lT9| , pO| , which argues that many nuclear 
properties are controlled by low-lowing moments of the Hamiltonian. Given a set of one- 
and two-body interaction matrix elements (the same input for a shell-model diagonalization 
or path-integral calculation) one can readily find in the literature straightforward formulae 
for first through fourth moments |T9|-pT|. One early result of statistical spectroscopy was 
that for a finite model space, the total level density tends to be Gaussian [^]. Of course 
this means the model ignores intruders at high excitation energy, but it also suggests that 



the level density even at low energy is governed by the width, or second moment, of the 



Hamiltonian. Grimes and collaborators included 4th moment [^,0 in order to improve 



the description of the level densities. Recently Zuker suggested a modified approach ^5 
He argued that a binomial distribution rather than Gaussian arises most naturally from 
combinatorial arguments. In the limit of an infinite number of states the binomial becomes 
a Gaussian; but for a finite model space the differences are nontrivial as we shall discuss. 

This paper discusses in detail the implementation of a novel treatment of level density 
descriptions using statistical models introduced recently by Johnson, Nabi and Ormand 
||26| . Here we compare these models against a number of nuclei whose level densities can be 
fully computed by diagonalization of a two-body Hamiltonian in the sd- or ^/-shells. We 
also compare against Monte Carlo path integration calculations for nuclides where direct 
diagonalization is not possible. We find that in most cases it is the fourth moment which 
distinguishes among the exact, Gaussian, and binomial level densities, although we find 
cases in which the 3rd moment is important. 

Section 2 presents the formalism of our calculations. Section 3 shows some comparison of 
the Gaussian and binomial level densities with the microscopic level densities. In section 4 
we present a curious finding, that there is a relation between the fourth moment of individual 
particle configurations and collectivity, which suggests a possible way to use statistical spec- 
troscopy to model the low-lying, non-statistical collective structure of the nuclides. Here we 
partition the model space into subspaces, use the binomial distribution to construct partial 
level densities and finally sum them to get the total level density. These improved models 
depict a more realistic picture of the microscopic level densities. We conclude and summarize 
in Section 5. 

II. FORMALISM 

Nuclear statistical spectroscopy championed by French and collaborators, is built 

upon moments of the nuclear Hamiltonian H in a. finite many-body space. The first moment 



4 



or centroid is 

^^, = H^{H)■ (1) 

all other moments are central moments: 

lir.^{(H-HY). (2) 

As one computes higher and higher moments, one naturally regains more and more informa- 
tion about H . But higher moments are difficult to calculate. An alternative is to partition 
the model space into subspaces Sa, using projection operators Pa = J2aea W) {(^\- We call 

^in{a) = {Pa{H-Hy) (3) 

a configuration moment, and general formulae can be found in the literature for configuration 



moments up to n = 4 ||19| , ^ . The total central moments can be easily found by summing 



over the configuration moments. 

The density of states can be formally written as 

p{E) =tT6 [e - h) . (4) 

One can also introduce configuration densities (sometimes called partial densities) 

Pa{E) = ti P J [e~h), (5) 

to which we will return later. If the many-body matrix elements of H are randomly dis- 
tributed (specifically, if they belong to a Gaussian Orthogonal Ensemble) then the density 
of states of H follow a semi-circle distribution. If H however is a two-body operator then 
the many-body matrix elements are correlated and the density of states is Gaussian, or at 
least nearly so . 



The Gaussian and semi-circle distributions differ in their higher moments. Let the scaled 
moments be defined by 

m„ = /i„/(/i2)"/^ (n>2), (6) 



as the width jj2 provides a natural energy scale. The scaled fourth moment is 3 for a 
Gaussian but is 2 for the semi-circle distribution. A Gaussian is reasonable good starting 
point as it is already close to the actual level density. A common generalization is to expand 
in a Gram-Charlier series using Hermite polynomials this unfortunately can lead to 
negative level densities. Another generalization is to extend the Gaussian to a function of 
the form exp{—aE^ — PE^ — ■yE'^ . . .) p7| , p4| but the relation between parameters 7 
and the moments is not amenable to a simple analytic formula. 

Recently Zuker gave a combinatorial argument that one should use a binomial distri- 



bution rather than Gaussian to approximate level densities We only summarize his 

approach here. Consider the binomial expansion 

(1 + A)^ = EA^(^) (7) 

Now interpret this binomial expansion as a distribution of levels, specifically, that A*"' (^^^ is 
the number of states at excitation energy E^ = ek, e being an overall energy scale. Because 
we can write (^^j with gamma functions, one can easily generalize to the continuous case, 

n(F ) = A^^/^ ^{Emax/^ + 1) 

^' nEje + l)TiiE^a.-E,)/e + l) 

where E^ax = ^N. Although we began with as an integer, it no longer has to be. 

The binomial distribution has several advantages. In the limit A = 1,A^ ^ oo one 
regains the Gaussian. For A 7^ 1 the distribution is asymmetric, because of a non-zero third 
moment; Zuker makes a combinatorial argument that the third moment can be as significant 
as the second moment. Finally, unlike most generalizations to the Gaussian, such as adding 
a term to the argument of the exponential, one can easily compute the moments of the 
binomial. The total number of levels, which is in effect the 'zeroth' moment, is 

d=il + \f, (9) 

while the centroid and width are given by 
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and the scaled third and fourth moments are 



Ref |25] does not give the fourth central moment. For Gaussians (A^ — * oo) 1714 = 3. For 
most binomials and for shell model diagonalization, the scaled fourth moment is less than 
3, a typical value being around 2.8. (Note that the above moments are exact for discrete 
distributions but are only approximate for the continuous distributions, due to errors in 
integration. For large A^, however, they are very good approximations.) 
Using Stirling's approximation, and a few others, Zuker arrives at 

piE.) - ^/^exp (-(AT - 1) In ^ + In ~^- ) + In a) . 

(14) 

The key parameters of the binomial are the order and the asymmetry parameter A. 
If A = 1 then the binomial is symmetric and has ms = 0; if A 7^ 1 then the binomial is 
asymmetric or skewed. The skewness can be positive (ms > 0,0 < A < 1), or negative 
(ms < 0, A > 1). Zuker suggested that the order of the binomial, A^, be fixed by the 
dimension of the model space. In that case A^ and A are fixed by solving eqns. (^ and (p!2D 
simultaneously. This we consider to be the 'standard' binomial. We note, however, that one 
could instead fix the order A^ by the fourth moment, and solve ([T^) and (|13D simultaneously 
instead, afterwards multiplying the entire binomial distribution by a constant so as to get 
the correct total number of levels. This we refer to as the fourth-moment scaled (FMS) 
binomial. After A^ and A are determined, the centroid and width simply fix the absolute 
scale of the distribution. 
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The central question of this paper is which distribution-Gaussian, binomial, FMS bino- 
mial, or some (other) improved model-best describes microscopic level densities? We use 
the low moments, in particular the third and fourth moment, to characterize our results. 



III. COMPARISON OF LEVEL DENSITIES WITHOUT PARTITIONS 

To test binomials as candidates for modeling level densities, we compare against exact 
shell model calculations. (All the densities shown in our figures are in fact state density, 
which includes 2J + 1 degeneracies. Strictly speaking, the level density ignores 2J + 1 
degeneracies, but the literature is often cavalier with this distinction). In this section we do 
not partition the model space into configuration subspaces, which will be considered in the 
next section. 

We considered full Ohu sd- and p/-shell calculations, and focused on several nuclides that 
could be completely diagonalized using the OXBASH shell-model code [EH]. In the sd-shell 



we used the Wildenthal USD interaction and performed exact shell model calculations 
for many nuclides covering the entire range. In the p/-shell we used the Brown-Richter 
interaction [^] and computed ^'^'^^'Yi. For still larger dimensions, we turned to Monte Carlo 



path integration calculations fT^HI^l- avoid the well-known sign problem |]T4| we fitted 



a schematic multipole-multipole interaction to the T = 1 matrix elements of the Brown- 
Richter interaction For all of these calculations we computed the exact spectroscopic 
moments, using the same single-particle energies and two-body interaction matrix elements 



were used by OXBASH. As far as we know this is the first time the exact formulas of |21 



for third and fourth moments have been computed for a general case . 

We ask: Is there a significant difference between exact, Gaussian, and binomial level 
densities, and is the difference driven by third or fourth moments? First the moments. 
Table 1 shows the scaled third and fourth moments for the exact distributions (the moments 



from exact shell model eigenvalues and those calculated using pT],|3T| should and do agree) as 



well as fourth moment for 'standard' binomials for some selected sd- and ]9/-shell nuclides. 
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(Remember that the third and fourth moments for a Gaussian are and 3, respectively. ) 
The third moment is generally small, except for the case of ^°F, ^^Ne, ^^Na and ^^Ar in the 
sci-shell and ^^Mn and ^'^Fe in the p/-shell; we found a strong correlation between nontrivial 
ms and 7^ 0. Because the standard binomial parameters and A are found by solving 
Eqn. (^) and (|T^ (that is, constrained by the total number of levels and the third moment 
), the fourth moment, m^, is not constrained. We see from Table 1 that on the average, 
the (standard) binomials tend to overestimate the scaled fourth moment, 7714. We further 
note that whenever the magnitude of the scaled third moment is > 0.05 (which we regard 
as a significant value), the difference in the exact and binomial scaled fourth moment is 
considerable, albeit less so for the p/-shell nuclides. 

Fig. 1 shows two sample sd- and two p/-shell nuclides. The histograms are the result 
of direct diagonalization. The figure shows the exact spectrum, the Gaussian and the two 
binomials (symmetric and asymmetric). For the symmetric case we forced the asymmetry 
parameter A to be 1 (the fourth moment did not change significantly). Because the symmet- 
ric and asymmetric binomial distributions are very similar for all these cases, we conclude 
that it is the fourth moment rather than the third moment that drives the difference. For 
^^Ti, in fact the symmetric and asymmetric binomials are indistinguishable. We remind the 
reader that our goal is a factor of 2 accuracy in modeling the level density. 

But what if the scaled third moment is significantly different from zero (i.e. if the 
absolute value is > 0.05)? Intuition dictates that the third moment should play a role for 
such cases. Figure 2 shows the case for nuclides with large asymmetry. Again we chose 
two sd- and two p/-shell nuclides. The open circles in Fig. 2 show the Monte Carlo shell 
model calculations along with the uncertainties. Here there is a clear difference between the 
symmetric and asymmetric binomials as expected. For these cases there also exist a large 
discrepancy between the binomial and exact m^. We therefore also plot the FMS (fourth- 
moment-scaled) binomial, which has the correct 777,4. The FMS binomial has clearly a better 
behavior especially for the p/-shell nuclides. Unfortunately the scaling also leads to a cut-off 
(unlike a Gaussian, the binomial has sharp cutoffs and infinite slope at = 0, Emax) which 
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is too high for the case of Ne. 

The prehminary lessons that we learn: binomial distribution do model shell-model level 
densities better than Gaussians, but in many cases it is the fourth moment rather than 
the third that leads to a discrepancy. Nuclides with T3 7^ have larger asymmetries, 
however, and must be accounted for. A scaled third moment of about ±0.05 is non-negligible. 
While an obvious solution is to adjust the binomial to reproduce the exact fourth moment, 
the binomial's abrupt cut-off in certain cases, leads us to generalize our approached to a 
partitioned model space, described in the next section. 

IV. SUM OF PARTITIONED BINOMIAL (SUPARB) LEVEL DENSITIES 

To improve our results we partition the model space into configuration subspaces. One 
expects partitioned models to better approximate the higher moments p2|j33[| (for example, 
Pluhar and Weidenmiiller had a sum of distorted semi-circles and got much better 
modeling than a single semi-circle). We also saw in the previous section that the FMS 
binomial gets the fouth moment correct but stops abruptly. Partitioning of the model space 
into subspaces solves this problem. 

We partitioned the model space by single-particle configurations, such as 
(131/2)^(0^3/2)^(0^5/2)^, for the very simple reason that analytic formulas exist for the as- 
sociated configuration moments. Other possibilities would be to partition by J, T or by 
SU(3); unfortunately no analytic formulas exist for these partitions. We did consider J, T 
partitions, computing the moments by hand (the OXBASH shell model code projects the 
Hamiltonian onto J, T partitions) but we did not find significant advantage. SU(3) would 
be an attractive partitioning , as it would naturally build in quadrupole collectivity, but 
even by hand this is difficult for a general two-body interaction. Experts should note that 
protons and neutrons orbits were considered distinct, which would allow us to easily project 
exact isospin. 

To compute the level densities, we take the following steps: (1) We compute the config- 
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uration moments up to 4th order. (2) We model the partial density for each configuration 
as a binomial (either standard or FMS). The binomial parameters and Aq,, as well as 
the overall energy scale and centroid, are fitted to the configuration moments for the a*'* 
subspace. The binomial then gives a partial density. (3) The partial densities are summed to 
yield the total level densities. Because of these ingredients, we will refer to our approach as 
SUPARB (SUm of PARtitioned Binomials) level densities. If we model the partial densities 
for the configurations as a FMS binomial, i.e., if we fix the order Na by the fourth moment, 
and solve ([l^) and ([13|) simultaneously for each configuration, we call the resulting sum as 
SUPARB-FMS. 

Fig. 3 models the state densities for four S(i-shell nuclides, ^^Ne, ^'^Mg, ^^P, and ^^S. Here 
we plot the SUPARB and SUPARB-FMS state densities against the exact ones. We see 
that not only the SUPARB-FMS densities model the exact state densities some what better, 
some structure in the low-lying states is also revealed. Such structure cannot be described 
by simple models such as either a Gaussian or a binomial. The structure is particularly 
very profound for the low-lying states of ^^Mg and ^^S. There were however cases where no 
appreciable difference between the SUPARB and SUPARB-FMS state densities was found. 
Fig. 4 shows four such cases: ^^Ne, ^^'^^Na, and ^^Mg. We see that the comparison is again 
very good (within a factor of 2) and that there is no appreciable difference between the 
SUPARB and SUPARB-FMS state densities. 

The computation time of moments, scales as (number of partitions) (number of j- 
orbits)^". That is, calculation of fourth moments requires much more CPU time as compared 
to the third one. SUPARB-FMS is sometime a clearly superior model but there are also 
cases where SUPARB alone gives equally good results. Two important questions then arise: 
(1) Is it possible to determine the cases where SUPARB state densities are sufficient, merely 
by looking at the configuration third moments in order to save time? and (2) What are 
the reasons behind when SUPARB-FMS gives better results? When do we see the low-lying 
structure in SUPARB-FMS? 

The answers to these questions do not seem to be simple. There are clues, however, to 
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the choice between the two SUPARB densities. One clues hes in the spread of the scaled 
third configuration moments. If there is a spread of these moments both above 0.05 and 
below -0.05, SUPARB-FMS is likely to produce better results. Fig. 5 shows the spreads in 
the scaled third configuration moments ma (a), where a labels the configuration, for some 
selected nuclides from Figs. 3 and 4. The ordinate is the associated configuration centroid 

which indicates the average excitation energy of states in the subspace a. (We have 
neglected a few of the configuration third moments in all the four cases with values much 
less than -0.05.) The bulk of the m3(a) lie in the range —0.2 < ma < 0.2. Any m3(a) which 
falls outside the rectangle is to be taken as significant. What we see is that for the graphs 
in the left column there is a spread of ms around the rectangle. For the right column we 
do not see any m3 > 0.05. And we see that for these two cases (^^Ne and ^^Na), standard 
SUPARB densities do equally well (see Fig. 3). So for these cases there is no special need 
for computation of the scaled fourth moments (similar is the case for ^^Mg and ^^Na). 

There can be a problem for high positive values of ms, A — > 0, and the binomial m3 — »• 
-7=. So, for example, if the dimension d = 500, by Equation (|T2]) the maximum allowed 



by the standard binomial is 0.16. For any value of exact > 0.16 the standard binomial 
fails. Frequently we found large values of ms for configurations even of small dimension. 
Fortunately, for fourth-moment-scaled (FMS) binomials, this is not as limiting (we fix the 
order by the fourth moment, solve ([l2|) and (|13D simultaneously, and afterwards multiply 
the entire binomial distribution by a constant so as to get the correct total number of levels) 
and the technique can be applied to all cases. 

Regarding the answer to the second question, the reply lies partly in what we already 
described previously. Empirically we find that a spread of exact m^la) around ±0.05, is 
likely to reveal some structure. Again if we look at Fig. 5 we do see that for ^^Ne and '^^P 
there is a spread of above and below the rectangle; simultaneously. Fig. 4 shows nontrivial 
structure at very low excitation energy for these nuclides. The same pattern occurs for ^^Mg 
and ^^S, suggesting that the two events are correlated. 

Yet another clue lies in the m4{a) shown in Fig. 6. We note that for cases where SUPARB- 
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FMS densities reveal some structure, the standard binomials systematically overestimate the 
configuration scaled fourth moments. One clearly sees that for ^^Mg and ^^S, the (standard) 
binomial 777.4(0;) are biased high compared to their exact values. On the other hand, for 
^^Na and ^^Mg, relatively structureless at low excitation energy, the binomial 7774 (a) 's are 
not biased relative to their exact values but fall both high and low symmetrically. 

Let us summarize what we have said so far for these improved statistical models of level 
densities. SUPARB-FMS density can be a better choice, not only in modehng the secular 
behavior (within a factor of two) but also reveals some structure in certain cases, specially 
for the interesting low-lying region. It however, requires a larger CPU time which can be a 
big problem when modehng thousands of cases. SUPARB density is still better compared 
to a single binomial density but for large, positive 7773 (which is very common for partitioned 
subspaces) they are constrained by the dimension of the configurations. 

V. CONCLUSIONS 

We have outlined a theoretical approach to level densities that is both microscopic in 
origin and also computationally tractable. We analyzed our results in terms of the third and 
fourth moments. For small asymmetries, the fourth moment dominates. For cases where the 
scaled third moment is significantly different from zero (with magnitude greater than 0.05) 
the binomials are characterized by both the third and fourth central moments. The FMS 
binomial gives the best description of the exact level densities, specially when the nuclides 
have a large scaled third moment. There is an improvement over the asymmetric standard 
binomial by as much as a factor of 2 in the modeling of the low-lying exact level densities. 
Our study also shows that the Gaussians continue to give a good estimate of the exact 
spectrum for the level densities (within a factor of 5) and this comes quite handy since they 
just require the centroid and width of the Hamiltonian as input parameters. They might be 
used for an initial estimate of level density and then the interesting cases can be followed 
by still refined models presented here. 
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Partitioning the space considerably improves the binomial approximation to the level 
densities. The SUPARB-FMS has some in-built structure and follows the trend of the exact 
level density specially in the collective regime. The degree of precision of the present ap- 
proach will give astrophysical nucleosynthesis calculations a much better predictive power. 
Future application to higher shells will be hampered as much by our ignorance of the effec- 
tive two-body interaction as anything else. SUPARB can also be easily generalized to the 
calculation of spin-cutoff factors, and we are planning to make application to nuclear heat 
capacities and estimate of contamination by spurious center-of-mass motion. 

This work was performed under the auspices of the Louisiana Board of Regents, contract 
number LEQSF(1999-02)-RD-A-06; and under the auspices of the U.S. Department of En- 
ergy through the University of California, Lawrence Livermore National Laboratory, under 
contract No. W-7405-Eng-48. 
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Table (1): Scaled third and fourth moments. The exact fourth moments are compared 
with those of standard binomials. 
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FIG. 1. Comparison of exact shell model state densities (histograms) against Gaussians and binomials 
for sd- and p/-shell nuclides. 
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FIG. 2. Comparison of exact shell model state densities (histograms) against Gaussians and binomials 
for sd- and p/-shell nuclides. The circles represent the Monte Carlo path integration calculation in a full 
Qhio basis along with the uncertainties. FMS = "fourth moment scaled" (see text) . 
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FIG. 3. Comparison of exact shell model state densities (histograms) against sum of partitioned 
binomial densities (SUPARB) and sum of partitioned binomial densities which are fourth moment scaled 
(SUPARB-FMS). Note the structure in the SUPARB-FMS densities. 
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FIG. 4. Same as Fig. 3 but for four other sd-shell nuclides. Here there is little distinction between 



standard and fourth-moment-scaled binomials. 
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FIG. 5. Spread of configuration scaled third moments m3(a) as a function of configuration centroids 

h{a) (relative to the ground state energy). Any value of ms^a) outside the dotted rectangle represents a 
significantly high value. 
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FIG. 6. Comparison of the low- lying exact (configuration) scaled fourth moments m4(a) with those 
calculated from the standard binomial. The dotted line is to guide the eye. 
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